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ABSTRACT 



A fundamental, non-linear, ordinary differential equation is 
derived for the flexural vibrations of a non-uniform beam, including 
shear . 

This equation is solved with the specializing assumptions of a 
constant shear coefficient and that the moment of inertia is directly 
proportional to the cross sectional area. 

A perturbation method is used and the eigen frequencies and mode 
shapes obtained are functions of a perturbation parameter e. These 
solutions are accurate up to the second order of e . 
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SYMBOLS AND UNITS 



ENGLISH LETTER SYMBOLS 
A Beam cross sectional area, I?* 

Lengthy combinations of other quantities (see equations (14) 
and (14a) 

2 

C-^ Proportionality constant, I = C-^A, L 

, -2 
E Young s modulus, FL 

-2 

G Shear modulus , FL 

4 

I Moment of inertia, L 

k Shear coefficient, dimensionless 

L Beam length, L 

M Bending moment, FL 




s 


kAG(y 1 -ty) > F 


t 


Time, T 


V 


Shearing force , F 


X 


Distance along beam, L 


y 


Beam deflection (positive downward), 


GREEK 


LETTER SYMBOLS 


ex 


, radians 


© 


L 



L 



r 

<T 



Mass density, FT L 

Tefe, F " 1 



f 

U) 



Slope of beam centerline due to bending 
Frequency of vibration, T ^ 
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OTHER SYMBOLS 



/ 





Symbols not seen here may be found in Appendix A. 
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INTRODUCTION 



The Timoshenko beam equations have been well established [5]*. 
Neglecting rotary inertia, these equations become: 



The left hand sides of both equations (1) and (2) are functions of 
y and i|r . Because of this common dependence, these equations are 
referred to as coupled. Uncoupling of equations (1) and (2) involves 
their reduction to one equation in which either y or \|/ , but not both, 
appear . 

It has been shown that the correction for shear is much greater 
than that for rotary inertia [10]. The effect of shear becomes 
increasingly important as the frequency of vibration is increased 
[9]. Therefore, the effects of shear are worth investigating. 

Figure (1) is a free body diagram of a beam element. Note that 
the quantity i Q Figure (1) is the slope of the beam 




( 1 ) 



( 2 ) 



centerline due to shear deformation. 



Figure 1 



X 



M 




*List of references appears on page 53 
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The shear term in equation (1) contains the shear coefficient, k. 
Timoshenko defined the shear coefficient as average shear stress over 
maximum shear stress. This shear coefficient, then, is a function of 

the cross sectional area of the beam. Timoshenko found it to be con- 

\ 

stant for general classes of geometrical shapes. However, there have 
been other studies investigating this shear coefficient which have 
disagreed with Timoshenko's values. In 1966, G. R. Cowper [4] 
derived expressions for the shear coefficient which were functions 
of the geometry, Poisson's ratio, and the dimensions of the cross 
section in some instances. R. D. Mindlin and H. Deresiewicz [9] 
maintain that Timoshenko's shear coefficient is a function also of 
the frequency of vibration. The following table contains some 
representative values for the shear coefficient. These values 
were determined by different methods and different assumptions. 
Therefore their comparison is difficult. 

Table 1 (From Cowper [4]) 





k for rectangle 


k for circle 


Timoshenko 


0.667 


0.750 


Mind lin 


0.822 


0.847 


Cowper u = 0 


0.833 


0.857 


u = .3 


0.850 


0.886 


c 

ii 

Cn 


0.870 


0.900 



There seems to be a conflict of opinions in the literature on the 
shear coefficient. This paper develops a general equation with k 
as a function of beam length but all solutions treat k as a con- 
stant. This is justified by the fact that for the class of areas 
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investigated, that is, where the moment of inertia is directly 
proportional to the cross sectional area, there is general agreement 
that k is a constant. 

The remainder of this paper will be devoted to the uncoupling of 
equations (1) and (2) and the subsequent solving of the uncoupled 
equation for its eigen frequencies and associated mode shapes. This 
uncoupled equation will be general with respect to cross sectional 
area, moment of inertia, and shear coefficient. That is, these 
quantities will be functions of the distance along the beam. 

Solutions of this type would be of use in optimization of design. 

It is often desirable to minimize deflection, bending stress or modify 
frequency of vibration in beam design. By allowing the cross 
sectional area of the beam to vary with x it may be possible to 
reduce deflection and bending stress or to modify the frequency of 
vibration through appropriate selection of an area function. In 
this thesis the area function will be defined in such a way that 
only one dimension of the cross sectional area be a function of x. 

For the sake of comparison to existing solutions, this variable 
dimension must vary in such a way as to keep the beam volume con- 
stant. That is, the beam volume must be the same if the variable 
dimension is a constant or a function of x. This stipulation will 
become evident in the derivation of the area function. 
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DEVELOPMENT OF THE GOVERNING EQUATION 



With the substitutions : 

s=kAG(y-« 



equations (1) and (2) become: 

[El 'f , ] / +S= o 
f Ay = s' 

where * - ^ and * = 

Solving equation (5) for ijf : 

f- y- scr 
r= y -(s<r)' 

Lp- 1 - y" / -fe<r) // 



(5) 

( 6 ) 

(7) 

( 8 ) 



(9) 

( 10 ) 

(ID 



Equation (7 ) is : 

ElfVETV+s = o (i 2 ) 

Substitution of equation (10) and (11) into (12) gives : 

El[v"'-s"<r-^s'<r^]+EI / Q/"-s / <rJ-vs(l-(ei<T')j=0 <«> 

Substituting equation (8) into (13) : 

EI[/'- T Wi- 2<r f , ^]- f El'(^"- a f'AS]+s[i-(EI<r')j= o 04) 

Writing equation (14) as : 

A* + sB*=o ( 14a > 



thus , 

S= (15) 

where BHo (see Appendix B) 
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Differentiating equation (15) with respect to x gives: 

rAj- 

fAj>B* 1 +A*'B*-A‘B* = 0 <l6) 

This form of the equation does not seem to have been previously 
reported. It is clear that this single equation is fundamental to 
the study of vibrations of non-uniform beams, including shear. 

For the case of a prismatic beam equation (16) reduces to: 

EI&- + 

From this point on, the analysis in this thesis will focus upon and 
be confined to the case of simply supported beams for which the end 
conditions are : 

ftu = ft w « ft'Oxxo = (y * W = ° (»> 

These conditions prescribe zero deflection and zero moment at both 
ends of the beam. Solutions to equation (16a) are: 

g^smC^t+^j l n» ci8) 

g n (v)= sm A n x ' A n = < 19 > 

where : 

/ EX X* 

are the eigen frequencies of vibration for the uniform beam. The 
associated mode shapes are the corresponding g ! s. 

Equation (16) is completely general and it is clear that difficulty 
will attend attempts to find useful solutions which also satisfy 
equation (17). Therefore, in order to make the equation more tract- 
able, a specialization is made at this point. The area function is 
limited to widthwise variations so that the moment of inertia is 






(16a) 



13 



directly proportional to the cross sectional area and the shape of 
the cross section is limited to those for which the shear coefficient 
is constant. 

After expansion and combination, equation (16) becomes: 



{ - ~ a* 1 AA "? +QeCiyD - nf [©p<^j 

,, +se i + 

f /1A AiA'A w /.d/vT 



, -J 

- - 



(20) 



& [ee’c,y"-f3© i r «,^Q - 1 
t'fe^pooSi'J t £ f? a A' A^-y l ~A"'-6A . 'f ^ jg%rc,y 

A AA j [ 

+ A"&c,^ ''+ © r u> a y] - $£{3Q}uj^- 1 w>A " A j 

(3© ^a?-^ -t Oj'o^w/yf + 

[ = o 

Equation (20) is the governing differential equation for the 
vibrating Timoshenko beam, neglecting rotary inertia, when the 
shear coefficient is constant and the moment of inertia is directly 



proportional to the cross sectional area. 
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SOLUTION OF AREA FUNCTION 



Because of the high degree of non-linearity of equation (20), it 
does not lend itself to closed form solution. Accordingly, an 
approximate method of solution is employed. It is desirable to 
obtain a solution which permits verification for a prismatic beam. 
With these ends in mind, a perturbation expansion of the symmetric 
mode shapes is considered. 

y = sin tls + £ sin + (21) 

Physical arguments suggest that the substitution of equation (21) 
into (20) will generate a symmetric area function. Substitution of 
equation (21) into (20) yields a differential equation with the 
parameter e in it, and we will take e to be small. That is, a 
perturbation method of solution will be employed. This solution 
will be a perturbation from the fundamental mode of the prismatic 
beam solution. If the variable area function is appropriately 
defined, the perturbation method will provide a link to the 
prismatic beam solution. This link will be effected by setting 
the perturbation parameter equal to zero in the expression for 
the variable area. Let: 



A- Ao*^ € A| ■+• 6 a Aa 


(22) 


u>2 -he 


(23) 



In equation (22) A^ is presumed to be constant since the solution 
we are seeking is a perturbation with respect to the solution for 
a prismatic beam. 
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It is expected that the solutions of equation (20) will be 
valid only for small values of e since they will be perturbations 



from the prismatic beam solutions. 

When the perturbation expansions (21), (22) and (23) are sub- 

2 

stituted into equation (20), area terms of higher order than e 
are generated. These terms may be neglected because the perturba- 
tion, in this case, considers only terms of second or lower order. 
For example, consider the third term of equation (20): 

A' a A"= [eA' ^Ai] 1 [e A?+ e’Af) 

-+-J2 A^Tj [g 

3 

The lowest order of e generated here is e , therefore, all terms are 

2 2 

of higher order than p and the term A 1 A* may be neglected. 

When equations (21), (22) and (23) are substituted into equation 
(20) and common powers of e are equated, three linear differential 
equations in A^ , A^, and are generated. The equation for A^ is 
simply the constant cross sectional area case. 

E ■+’ © A 0 - f A 0 ^o y 0 = O 



For a simply supported beam: 



2 



CD 



0 




SC, PE ft 



is the leading frequency term in the perturbation expansion (23). 
The first order of e generates a linear differential equation 



in A^(x) . 





- a " [3 -©* 



-rc,*'- + a; QECtt'j* A, [2K, v ® (24) 

-3.fUJo + “ A c [EC,^^- \/, 

"f + ®(° >" + © f ^ 

This is a linear differential equation with known variable coefficients. 
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The only restriction on the total area function is that the volume 



of the beam must remain constant. Therefore: 

^Adxs Volume e^A.dx +6 a [ A A a dx (25) 

The perturbation parameter e is arbitrary thus defining a class of 
areas. Therefore, to keep the beam volume constant for all e , it 
is sufficient to require that: 

5 A A ' dx= L a ^ x=o (26) 

Assuming that A^(x) is piecewise regular and symmetric it may 
be represented by a Fourier series of the form: 

op 

+ Mlcos 

u- o U ' 



Equation (26) restricts the arguments of the Fourier expansion to 
even powers of . 

£ A, dx - £ S 0 [N + Me cos tT^jdx^O 

= ^-nrNi,(cos Ctt- I) + -fe* MtSincir] = O 



L=0 



Sin iff =0 ’ i- O, Iji, 3,- • 

COS ITT - I * o ; L = 0,a, — 

Therefore , for J A|ctx = O > it is sufficient that the arguments 
of the sine and the cosine functions in the Fourier series be even. 
Substituting : 

sm nx 

Vf,= sin 33g. 5 ©4 s igs. 

into equation (11) gives : 

a!', 1 K, cosoc + A[ # K a Si r\* - A,' k 3 cos»c - A, k 4 SI A Jk 5 si r»3« + ^sm^jc 27 ) 

Examination of equation (27) shows that only the cosine terms 
in the Fourier expansion remain in the solution for A-^(x) . if the 
sine terms were allowed to remain they would generate products of 
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cos a and cos na and products of sin a and sin na. These products 
can only be represented as cosine functions. The right-hand side 
of equation (27) contains only sine terms. Therefore the coefficients 



of the sine terms in the Fourier expansion must all be zero. On the 
other hand the cosine terms in the Fourier expansion generate only 
products of sine and cosine which may be represented entirely as sine 
functions through trigonometric identities. Therefore: 

A.« = Mn cos not 
r\~ o 

Substitution of this expansion into equation (27) gives: 

OQ 

K Sm(n°<+°<)+isin(rv*-*)] - (<* ^ smtfwh*] 

fcf j^sin(m+oO sm 3* 4 

There are only two values of n which can satisfy this equation: 

r\= zz 

n is taken equal to + 2 because, 

CoS 2t.o <. - cos C-*0 

Substituting n = 2 into equation (28) gives: 

M,Z,sir\3<* +- M,Z a sinoc « A 0 jk s s«n 3* + K 4 sm<*] 

Expansion of the terms K^and show them to be identically zero 
(see Appendix A). Therefore, equating coefficients of common 
arguments of the sine function gives : 




A, 00 = M,cos (29) 
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The unknown frequency is contained in the term Kg. 

M,Z,«A .k t =0 

- A of =uj, i [i+e(s\ ;! ]=o 

Therefore : 

cjf- O (30) 

Thus, the first eigen frequency has no contribution from the first 
order e solution. 

The second order terms of e in equation (20) produce a linear 
differential equation in A 2 (x) with known variable coefficients. 

AaK,COS<* + AlKa,Smo<-Ail<3COS«: -A* ^5,113* = 
sik\o< •+• si n 3°< 4- ^ 2 A SoC (3i) 

Equation (31) is similar to equation (27). To maintain constant 
volume, it has been required that 

J 0 L A i Wdx=O 

The solution of A^(x) , therefore, will be of the form: 

A*0O * 2. Mr\ cos nu (32 ) 

rv-o 

Substitution into the differential equation (31) gives 
The value of n = + 2, +4 are taken because cos x = cos (-x). 

A a 0O= M^COS^oC + MsCOsMoC (33) 

2 

This is the part of the area function whose coefficient is e . 
Substitution of equation (33) into equation (31) gives: 

MaZaSin* + [Kz,+ M 3 Z 3 ]sm 3oc sm Secs 

3oc + ^ 3 sin Soc (34) 

Z 2 is identically zero (see Appendix A), therefore: 
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The frequency o)^ is contained in the term £ . 

a * uj ^ [p +© f (j-rj - o 



^ ’ AoLp+©f(E)ii 

To find the unknown amplitudes, coefficients of sin 5a and 
are equated. 



( 35 ) 



sin 3a 



^3 = fyz* 

H a = ^ ~ KsZ3 

At this point the class of areas in terms of the parameter e is 

2 

determined to the second order of e. I s the contribution to 

2 

the first eigen frequency from the e solution. 

A = Ao+6 K Co s2^ + £*[MaC 0 ^* + ^ 3 COsH*] < 36 ) 



The frequency of the first mode of vibration has also been deter- 
mined to the second order of e . 

(37) 



rD+©eE) l 3 



Equations (27) and (31) are third order ordinary differential 
equations. They both have three boundary conditions imposed on 
them. For equation (27) the boundary conditions are: 



x=o * A,= M, 

x=l ; a,= M, 

^A,dx = 0 

The conditions on A ^ are: 

x. - o • A x = Ha-'-M 3 
x= L ; A^=M 3k + H 5 

( Aa dx = O 
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The solutions obtained for A^ and A 2 satisfy both the differential 
equation and the boundary conditions in both cases. These solutions 
are then, unique [13], 

The class of areas obtained here are symmetric about the midpoint 
of the beam for all values of e. For e = 0 it may be seen that the 
solution for the prismatic beam is obtained. 



21 



SOLUTION OF EIGEN FREQUENCIES AND MODE SHAPES 



Extensive algebraic manipulations have been carried out through- 
out the following work, For the sake of clarity of presentation, many 
long algebraic expressions are represented by symbols. The expansion 
of these symbols may be found in Appendix A. 

Using the results obtained for A^ and A^ , it is now possible to 
obtain equation (20) in the form: 

= o < 38 > 

For the simply supported beam, the boundary conditions are: 

; a = ° >3"=° 

* = L i " = o 

The procedure in solving equation (38) will be similar to the 

procedure used in obtaining the variable area functions A^ and A^ . 

That is, three separate differential equations in y A , y _ , and 

nO nl 

y will be obtained by equating like powers of p , on either side 
n2 

of the equation. Solutions to these three equations will give all 
the eigen frequencies and deflection modes in terms of x and the 
perturbation parameter e . 

The notation used will be as follows : 

= + (39) 

C On = C 0*0 + €GOjt, + &U)n3. (4o) 

The first of the double subscripts refers to the eigen frequency 
or deflection mode. The second subscript refers to the order of e. 
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The zero order of 0 gives the constant cross section equation. 



+& - fW 0 sy no = O (A1) 

&«o)x»o * Cj^) k . u - o 

The mode shapes are of the form: 




The amplitude , B^, may be determined from initial time conditions. 
Substitution of equation (42) into equation (41) gives: 

BnoSm laX[EC, In -Qf^o A* ~ f = O 



Therefore : 

OJno~ 

This result, 



f li + © wj 

for n = 1, checks with 



(43) 



the earlier result in Section 



IV. 



Equation (24) may now be expanded into a differential equation 

in y . By equating the coefficient of e in equation (38) to zero, 
nl 

we get: 

+Q a ^',' -Q 3 ^> = Qi+ sm a* +Qs-sia3ot c-os n«. 
- cos 3 -ol 51 a rv«?c 

Applying trigonometric identities : 

+ Qa. \j nl - Q 3 \f n , ” Qq 5/ A ft* (44) 

“ a [jQs *+0^ s<ft( ^ £c?,y ~ nO'W'2) oc 

where , 

— T5no P 0 XryJ 

If were non-zero, its contribution to the solution of equation 
(29) would be non-periodic [7]. Therefore, its value must be 
chosen such that it vanishes for all values of n. 
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A Fourier series is assumed for y _ . The coefficients of all the 

nl 

cosine terms in this series must be zero because the right hand side 
of equation (44) contains only sine terms. After substitution of the 
Fourier sine series into equation (29), only three sine terms have 
non-zero coefficients. 

\j M = Bm sm(n.-2)o4-hC n ,sinaot +Dm sinOua)* (45) 

The requirement of - 0 for all n makes C ^ = 0 for all n. 

Substitution of the assumed shape for y n into the differential 

nl 

equation (44) gives : 

+• DrtijQ^-Qa^-QjJ 

Sm(n^o< =Q 4 sin°< +^ y -Q 4 Jsm(n+;k-^jQ s i-Q t JsinCri-l)<?c 

where , 

P= (A-A)JT- 

<^= (Af^) 

When n = 1, the sin(n-2)a term becomes -sin a. 

Therefore : 

for n = 1 shows that it is 
Therefore : 

COn =0 

2 

This confirms the earlier finding in the A^(x) solution of = 0. 



[Qs -%-%] = O 

Expansion of the quantity ijQs+Qt 
identically zero (see Appendix A). 

f; 

q m =o 
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To maintain a periodic solution, must equal zero for all 



other values of n. Therefore: 



uJ A ,-o } l 



The amplitudes B^ and D are found by equating terms with 
the same sine function arguments. For n = 1: 

B„= CQ-t - 5P - 

Qt p^-^p 2 -- Q* 



Dn = 



Qs ~ Ql 



^LQ,a‘ f -Q a a z -Q<3 

The expansion of shows it to be identically equal to unity. 

This checks the earlier substitution of: 

+ +< ^x > y,- , a= I 

For other n: 

B,v, = jlPs+Qj] ) A > I 






Table 2 



n 


Amplitudes 


Frequency 






T^ru 


Ca-J/ii 


/ 


o 


1 


o 


n>l 


- CQs+Qd 


(o^-Q«n 


o 


Jlap^-o^p’-Q^l 





cu _ is the contribution of the first order e solution to the overall 
nl 

, 2 
eigen frequency . 

U-)r\ ~ LUno + & + S ~ Oj? A ^ 

B and D . are the amplitudes of the sine functions which 
nl nl v 

contribute to the overall equation for the mode shapes. 

^0+ £y>\i 
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where , 



The second order e 

Ao[Q.bi“+Qo- 



equation in the unknown y is : 

n2 

c^i sm(n4aW + 6a3mCa-^oc 



+ <£3 31 nC n+q ) oc 4- 614 2>iaC n~H)<x 4 <$5 strv. a.<* 



( 46 ) 



The assumed shape of is therefore: 

^)ax z B rva SiaCa*v5.')cx*+T)»\3L siaCa _ ^)c< •+ £rva sm(r\+if)<* 

+ Fr\iS/aCa-'4)cx r G?n 2 SIA aoc 

The term 6^ is : 

&5- BrVQ f Ql+0 XnJ &<o 

As was the case with the first order c solution, this term may not 
be allowed to persist if the solution is to be periodic. This 
condition defines the frequency and demands that G ^ = 0* 

For n = 1 , the §2 sin(n- 2 )a term becomes - §2 sina. 

Therefore : 

a _ <£3- ^6 

B,<.pD+©®)*jA» 

Similarly for n = 2 , the 6^ sin(n- 4 )a term becomes - 6^ sin 2 a 
and combines with the 6^ sin 2 a term. Therefore: 

uf~ - 

The amplitudes 5 D 2 5 ^ n 2 ’ anc * F n 2 are ^ oun( ^ 5 as in the fi rst 
order e solution, by equating coefficients of the sine functions 
with like arguments. Care must be taken, however, noting combina- 
tions for certain values of n. Substitution of the assumed form 

of y into the differential equation ( 46 ) gives : 
n2 
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Bvva, (0»^ - Qa <^ 2 ~ sinCnfa^ +IW [Qi p 4 -Qip i -G^ 3 J 

St^Crv -<*)<* ■+■ EnafQtg/' 4 * Qa^ /a -Q-J SirvCl\+M)o< 4- 

BvjQ, ^m(^K+(SaSlA(n-A)0C 
4-^3 f^qsmCn-Hicii-^ssta k <*]~ a ? 

where , 

p'- (h-4)l 

s (ft vH) ^ 

For n = 1 

The expansion of 6^ - 6^ shows it to be identically zero for n = 1 
(see Appendix A). The expression fog^-oaf- q 3 1 is non-zero. 
Therefore [B^ " F^lj w hich is the coefficient of sin 3a for n = 1, 
is zero. 

Q*(r£ ) x - Oz\ - - ia.) a* 

2 

From the earlier calculation of a) > we know that [6,- - § 2 ] =0. 
Therefore : 



D,z=0 



^ ,a lPi(^ M - Q»(®5 p ) ~ Q»J ~ 

Expansion of Q^-Qj"] At 



^/ a & 

shows it to be identically 



equal to one for n = 1 (see Appendix A). 




This checks the earlier substitution, the the A 2 solution, of: 

^ = 1 sm 5 tdl n»i 

M« 9 
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o 



n> l 



Expansion of the expression for oo^ a lso shows it to be consistent 

2 

for the value found in the solution. 

For n ~ 2 the sin(n-2)a term vanishes. Similarly, for n = 4 the 
sin(n-4)a term vanishes. Therefore: 

o 

Fhj.= o 

For all other n. 

n zsr — i i ft ">2. 

rvx ” [Q,p M - Ao 

P*-.- — , > A <5- I 

- Q^'*- Q-a3Ao 

c Sa — — , . rv**H r\k\ 

2 nx ~ D5.p /l -<iV a -Qi) Ao 

cd o is the contribution of the second order e solution to the overall 
nz 

2 

eigen frequency cd^ . 

— LOno + € -f- 6^na 

B 0 , D OJ E and F 0 , are the amplitudes of the sine functions 
nZ nZ nZ nZ 

which contribute to the overall equation for the mode shapes. 



where , 

>*= 3 na. S i n0i+2)o< -t-TV^si nCnriU +E na si a(n+M )* + Fvi a Si a (ri-M ) «c 

The mode shapes for the class of areas derived in Section IV are: 

\jK=B ft0 SlA!Vg< +e[Brv,Sllf\ p* +Dm sm^xj 4- 

^ {Bna. Sift cj* +DaiSia px +Eni sm^x +Fru5in p'xjj ^ 7 ) 
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( 48 ) 



The eigen fr equencies are : 

UJ*. 



_ ^*5in obC Xv\ ^ X>A 



Br\o fQl +• © inj" 



where , 



X, - (Sx-^)^i “a? 

Xi* -<s 



4 



n >cZ 



is an amplitude which depends on the initial time conditions 
of the problem. 

By setting e 5 0 in equations (47) and (48), the solution is 
reduced to the constant cross section solution. This is expected 
because when e = 0: 



A - A 0 = constant 



By setting n = 1 in equations (47) and (48), the values of 
and are : 

\|.= smiH+esm Sin Six*. 

v -'* L. L. 

These results are identical with the results and substitutions from 
Section IV, the solution for the area function. This checks the 
great amount of algebra which was performed in obtaining formulas 
for the eigen frequencies and mode shapes. 
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CONCLUSIONS 



A fundamental non-linear ordinary differential equation, equation 
(16), has been derived for the flexural vibrations of a non-uniform 
beam, including shear. 

Equation (20) is the particular form taken by equation (16) under 
the specializing assumptions that moment of inertia is directly pro- 
portional to cross sectional area and that the shear coefficient is 
constant . 

Tables 2 and 3 contain approximate solutions of equation (20) for 
all eigen frequencies and mode shapes in terms of a perturbation 
parameter e. The perturbation solutions are accurate up to the 
second order of e . 
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EXTENSIONS 



The preceding results could be applied to beam optimization design. 
Using the formulas of the preceding sections, the deflection may be 
obtained as : 



It may be desirable to minimize the deflection y. For example, 
minimization of the amplitude of vibration could cut down on noise 
level. Bending stress is a function of the deflection and this 
quantity could also be minimized. This minimization process could 
be approached through iteration on the perturbation parameter e . 

A value of 0 would be chosen and y(x,t,e) would be evaluated at 
increments over the beam length over appropriate time intervals. 

The maximum value obtained from this process would be saved. This 
process would then be repeated for other values of 0. The minimum 
value of the maximums noted for each iteration would give the optimum 
value for 0. 

It must be noted that the equations developed are valid only for 
small values of 0. 

Another possibility for further study would be the application 
of this method to the forced vibration case and subsequent optimization 
of design for deflection, stress, or frequency. 

An experimental investigation to compare with theoretical results 
would be useful. Problems would probably arise in the manufacture of 
of the beam itself and the development of devices sensitive and accurate 
enough to detect variations from the solution for a prismatic beam. 
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APPENDIX A 



EXPANSION OF SYMBOLS 

This Appendix contains the expansion of the symbols used in the 
main body of the thesis and those symbols used in Appendix D. Symbols 
equal to zero are shown. 

The expansion of the symbols used in the area solutions follows. 

K,-© 2 f^®-@EC.£f 

K*: 3& r ujf + e a l =c^ < , l ® 1 +ec,®f-eE'c,(5)‘ t 

i Qec rn 1 * • - © p <*>?<£ ) 2 J * o 

ecX 

K«- Spo -fou, a =o , =o) 

- ftoi - ©p 

K,= 

kv IE -puJo 1 - Bp uj} (S g-f 

KlO- ^ ~ o 

Kn- - ;?e a pw? -©©<:,<£ f 
ku= 9 -s©£Ci( 5)" 
k,v e>ujj®f-© ec,(EA+7©f^. a 
k«= -©EcfE) 1 -© 2 ?^* 

k lS « -3©pw, a - ©'pc^ffFf = o 

K, t = ©ec.w" - 3©p -©V^o^-eccs?;* 

k,r- ©*p^(¥) -©e-c,(a ^) 3 

k.s= ©>«->*(£) = o 
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z, = ^k,(c) 3 --Zk a (t) l +i<3(5l 
Z»= h k ,0 *2 k, £ jN- K,(t ) = o 
Z 3 = 3A k.CE^+SI^OE^ + ikj®) 
Z H = 33 k,<jf - 8 Kx®) 1 k 5 ffi-) 

(4 - Ao kj 

^a- A 0 ky - O 
Q>sr Aok, 

^rW.AoQc.o-GffK'J ^ N>= 

^s- Wi Ao(j ks - ft?)* k'O 

& = n, Ac (sMffi) ks +( 2 e L f k n3 

^7=N, J A 0 [(^fk, 2 - e-£) 3 k, -<3F)K{] 
&= N, 1 Ac[C= ! ?) i k, 3 +(3?) 4 ki-0 

^'o- N\ AcQ( ) ki^j = C5 
f I - i O’H' 3 ' "*■ @4 + ^V + 

f *= :£&& + @-* + <?j6. -QVa. + @yfe] 

^ [?<« ~ "^4— (3'^ - ^cQ 

Mi - S' ^ 

Mx= ^~ Z M3 ^ 5 
Ms' VZh 
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The following symbols were used in the solution for the eigen 
frequencies and mode shapes. 



Q» = £C, 

Os — PLOno 

Q h = BnoLf +©f ArvJOJ n/ -O ( C^Jni-O , /L- /1o£ 
Qs - Bko Ak&J 1 - ©^C, ^ 6?)*- j2 EC, l 3 J 

0&= Bko N 1 - ©ec, Af 4© a ^>co A o An 46TC, AnJ 

O"? — "Bno [jbC-t tit-peon* - Anj^O 

Q?*BnobeC, ft] 

C\^ =• Bnot^^fCUno “©EG AU" + © 2 (®U^ao Xn+BC, An] 

Qd “Bno An. - CuJo An3 

On — Bno C^4 ® f An J 

G?|-x= Bru Qf + ©f>p z ]u^ = O , p = 

Ora- Da» Cf+<3f 3*3 ca>* ~ o ^ q~ (n+2) *Er 
GJm - Bru DzC, p 4 - “ ©f 6 l >*0 p5 

O15 = Dai [jECiq 4 - fC/Jno ~ 0 f a^yfo C£ x ] 

On, — "Bno Cf > 4 ©C An.] kCV» / ~ (3 
On — Bru C=A EC.i p~] 

O* - Dai [^SC, 0_ 3 J 

Gm - Q* 

Qao-Bno Caer^Jo +©ec, ftl 

Q*i — Bno L ^ ©^U^no An- (p ©tCj An] 

Qaa.” Bno Ceerc* a^ - 0 ^* & +7e r mA<J 
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Cp 23 >- Br^ofeS’Ci +® 2 f>LUn l o] 

Q *4 = 8 m [©EC, - 3©f»<jUn*o -©V p 2 - EC, p aJ 

Q* s - Dm feec, a M - 3 ©fW c - © a fW° a 2 - <= c > at! 
Q *6 - BnoC^r^ * 4 A 2 ] = o 



Qsn 

Qas 

Q30 

■ 

<m 

4 z- 
♦* - 
%- 

4 s ~ 

4 ^ 

4 t = 

4 * - 
<K- 

4 «o- 

4 >n- 

4 >f^ 

^I 3 ~ 

4 m' 



= ~Cp«i 

- Bai Ct © f=D, p' -5 4 © 2 ^ LVno pn 

’ Dm C- ©EC, UJnb <^H 

= 1 A*!] ~ D 

• ~ C?iO 

- A d KaC®F) z Qql 

- AoMsDCtf ) 1 GJql 

- A 0 MaOcF) Qs^et)^Q><Q 
A 0 M 3 R^jQ? 4C t U") 3 QlO^J 
AoQu 
Ao c5(j 

Ao Q»3 - _ 

Ao Wi [AQih -Ctt) Q 34 I 
Ao Ni [pQis £n?) A Qas~I 
Ao /V 1 (S 1 ) 3 6?a $ - (^F) G?n jf 
Ao ^?isT[ 

Ao N*Llir) 3 Q3i -®)Qiq ~(§?) 3 Qa>3 

A1N^)*Q. o +(ZP)* Q>n_ 

A?W. a CW^Q« -(WO*. 
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Ai- i[l^+4>3 + 4>-T| 

A*- iC4i -^3 

A 3 =-i Qj> a - 4q - cj>cj + <f>„ + ^vi * 

i BK+4>q - A? - 4i° " + V* “ ^ J 

As = it C^As ”4? -4>q +4/o ~ + Ai3 +4 n^ 

Afe-*^Q"4& ■ + 4(0 -4** 4wj[ 

p'~ U-m) £ 

- tfi + 4) £■ 
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APPENDIX B 



THE CASE WHEN B* IS ZERO 

On page 10 an assumption that B* was not equal to zero was made 
Suppose B* did vanish: 

l -[eio -'] = o 

L^Icr'J '= I 

Integrating with respect to x: 

EI^=x+a 

If k = constant and I = C^A 

^ “kJq ^ 

~ j c , f dA r f X (/+ CA^dy 
KG, 3b A Jo 

'i%^ A = ;P + cxX + b 

A= exp[-^(f\ax4 lo)J 

This is not a symmetric function about the beam midpoint. Since we 
are restricting consideration to situations having symmetry about 
x = L/2, this case is excluded. 



38 



APPENDIX C 



NUMERICAL EXAMPLE 




Steel Beam E ^ 30x0* p S I 

Gi= 12x lO^ps i 

k*,£7 

r *WO(l> f/ft 3 



Ao= lm* 

Icbno- 1 1 A. - j A ^ I 

I=C|A ,C=ly,a= Jfem* 



A computer program was prepared that computes the first ten 

eigen frequencies and mode shapes for 0 - .05. The program also 

2 

computes and plots the variable area function cA^ + e A 2 for 
€ = .01, .02, .03, .04 and .05. The first ten eigen functions for 
e = .05 are also plotted. 

2 

The physical significance of 0 A^ + 0 A^ in this case would be: 

A- b(*)h= A q + £ At 4- A X 

Therefore, the width would be a function of x while the height, 
h, remains constant. 
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The first ten eigen frequencies are: 
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APPENDIX D 



INTERMEDIATE FORMS OF EQUATIONS 
Referring to page 25 the differential equation for A^(x) is: 

A 6 ( A* Kites* + A" Minot -AiK*cos*] = Ao 

^»> E - ' f +AoAf 

[ a ^"'1 + A, At{2EZ,tfl+ Ar^-je^^+GPC,^ 

+ g© s-c, y"| 4 A;*[-« z r w^'-ecc^f 
+7 ©|°^o A \/ 0 ^ - A, / A',"^©e'C,j o ''- & 2 pUJo'y^+h>A!"j6£0'^ 

- 3&pcu, 1 \j :} + ©^c^o \j 0 " 4©C, 

+ A, A,"f©©C,^?- 3© r 60o a ^+©V^V+ B^J-AoA," 
^ec, ^,'"4 © > - A , A,'" £©e:,y 'v © z fW y f 

Substitution of A-^(x), y^, y^, and y^ , and subsequent reduction 
gives : 

AoCA'a^^os^+AjkaSm^-AiK^cos^r]- ^s/a^+^siaB^ 

t^aSin 5c* -h^smetcos^* + (3s sin 3*x cos - 
s I ft 2>ot cos 3 oi + (2-? coso( cos<2o< s/n ^©t-f 

^5ift<*cos^c* + ^ 5(M sin 2 2 oi. - ^ cos^ s/n^oi 

After further reduction, using trigonometric identities, equation 
(31) is produced. 
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is similar to that for 



The basic differential equation for 

A 9 (x). The only difference is that y n 
L n2 

After substitution of A^ A 2 , y nQ3 and 
form : 




nl 3 



treated as the unknown, 
the equation takes the 



yn* -Q 3 y/jJ- 4 >iCo 5<5<^ sin t 4 > 2 CosHo< sinno/ 
+4*3 s,n cos n<x - cj >4 sm ^ cos n<x -k |>5 s/n nu m px 
+ <t> 7 Sinc^x -<j)% cosset smpx -4, cos 2c*: smcjx 
+4o sin^ cospx +4 m 5/n Pvoi cos<pc -f 4a cos^ sin2^ oo 5 /ix 
+4 g co5 2 Ax sm h(X + 4/11 3i n 2 Zo< si n nx 

Application of trigonometric identities to this equation produces 
equation (46). 
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